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Abstract 

Motivated by cosmic censorship in general relativity and string theory, we extend Christodoulou’s 
celebrated examples of naked singularity formation in the Einstein-massless scalar field system to 
include a positive or negative scalar potential of exponential forms, i.e., V{4>) = ±ex.p{2(f/K) 
with a parameter k. Under spherical symmetry and a continuously self-similar ansatz depending 
on K, we derive a 3-dimensional autonomous system of first-order ordinary differential equations, 
which incorporates the equations for massless scalar fields as a special case. Local behavior of the 
phase space is studied analytically with global solutions constructed numerically. Within the 3- 
dimensional solution manifold, we observe, for the negative potentials, naked singularity formation 
from nonsingular initial data for nf < 1. Meanwhile, transitions between solutions containing 
naked singularities and black holes are also identihed. However, when the potential is taken 
positive, numerical evolutions result in formation of black holes, but not naked singularities. 
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I. INTRODUCTION 


A (globally) naked singularity in general relativity is a singularity visible to distant ob¬ 
servers. The weak cosmic censorship conjecture ^ postulates that generically naked singu¬ 
larities do not occur from gravitational collapse of physically reasonable matter. In other 
words, singularities should always be preceded by a trapped region without future light 


cones expanding to inhnity. During the past few decades, on 


y limited progress has been 
4]. From 1980s to late 90s, 


made toward a general proof or disproof of the conjecture 
Christodoulou investigated global dynamics of the Einstein-scalar held equations with a 
vanishing potential in a series of papers. Remarkably, he established the hrst example of 
naked singularities that can evolve from regular asymptotically hat initial data in a suitable 


B- 


matter held [51 . Later on, he also proved that such occurrences are unstable and non-generic 
H. hence not violating weak coanac cenao.Mp, By conataucBon, Me faa.ily ot exact apace- 
times in [5[ is continuously self-similar and has two real parameters k [k in our notation) 
and a. Detailed analysis was given to the moment when the curvature singularity hrst ap¬ 
pears and the self-similar coordinate becomes singular as well. A naked singularity arises 
provided 0 < < 1/3. The asymptotically hat initial data sets are obtained by trun¬ 

cating the self-similar ones. Relevant to this work, Brady examined the same model using 


B 


numerical simulation [7i with more emphasis on critical behavior. Bedjaoui et ah further 
applied the self-similar construction to the Brans-Dicke theory Sj and showed that naked 
singularities can also emerge. In parallel to Christodoulou’s work on massless scalar helds, 
Choptuik discovered a diherent type of naked singularities, admitting discrete self-similarity 
instead, in critical phenomena through hue-tunings of initial conhgurations 


Q 


see also 

Since then, critical phenomena have been reported in many other systems (see the re¬ 
cent review llj), for instance, in the Einstein-scalar held system with inclusion of nontrivial 
self-interaction potentials (quadratic or double-well forms 12l-ll4l|). For more discussion on 


self-similarity in general relativity, see, e.g.. 




and references therein. 


In the context of anti-de Sitter/conformal held theory correspondence (AdS/CFT), there 
has been a series of discussions on possible violation of cosmic censorship in asymptotically 
AdS spacetimes ]^28|. Most attention has focused on the Einstein-scalar held system with 
a self-interaction scalar potential V (0). Unfortunately, for non-vanishing potentials (possibly 
negative, bounded or unbounded from below), the question of weak cosmic censorship in 
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spherical symmetry remains open 2^. Examples of naked singnlarity formation analogons 
to Christodonlon’s appear to be lacking. 

In this paper, we aim to show that Christodonlon’s self-similar model of naked singnlar- 
ities can be extended to inclnde exponential potentials. In cnrrent literatnre, exponential 
potentials (Lionville type) often arise from scale-invariant snpergravity t 


leories as well as 


the Kaluza-Klein dimensional redaction on curved extra c 
also been intensively studied in scalar held cosmology 27 


imensions 


24 


271 . They have 


28[. Most importantly for our 


purposes, the exponential forms are compatible with the assumption of continuous self¬ 
similarity, which makes our attempt viable in the hrst place. 

In presenting our model, we will follow Brady’s precedent of with certain modihcations 
(local compactihcation of the phase space) in order to exhibit the naked singularity solutions 
more prominently in numerical plots. When it comes to analysis of the equations and 
spacetimes, Christodoulou’s paradigm is respected in our treatment. In particular, we can 
embed his 2-dimensional solution manifold into our 3-dimensional manifold, allowing us to 
draw analogy along the way. 

The paper is structured as follows. In the next section, we introduce the self-gravitating 
scalar model and reduce the Einstein equations to an autonomous system of evolutionary 
equations. In Sec. Illll local studies are performed on critical and singular points of the 
phase space. Then in Sec. HVl we piece local pictures together and solve the equations by 
numerics in two parameter regimes k > 1 and 0 < k < 1, respectively. Afterwards, Sec. 
El is devoted to analysing the spacetime structure of the naked singularity solutions, and 
concluding remarks are made in Sec. IVll For interested readers in Christodoulou’s paper 
[^, the Appendix n contains our autonomous equations using his original notation. 


II. THE FIELD EQUATIONS 

We follow mostly Brady’s notation in for an easier comparison. To start with, the 
Lagrangian density is chosen as 


^ ^ (l “ - U(0)^ , 

V (0) = Vo exp , « 7^ 0, Uo = ±1 or 0, 


( 1 ) 

( 2 ) 
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where 0 is a real scalar field and V (0) its self-interaction potential with a parameter k. 
By a translation of 0, one can normalize the mnltiplicative factor Vq 7 ^ 0 to ±1. The 
energy-momentnm tensor is given by 


= V^ 0 Vi .0 


- (V„0V“0) + 1/(0) 


d/j,p 


( 3 ) 


The Einstein eqnations = 2T^j,) and the generalized Klein-Gordon eqnation are 


= 2 (V^0V,.0 V , 
V„V“0 = E'(0). 


( 4 ) 

( 5 ) 


Using retarded Bondi coordinates in spherical symmetry, we consider the following metric 
(g, g real-valned): 


ds^ = —g{u, r)g{u, r) — 2g{u, r) dndr -|- r^dhl^, 


( 6 ) 


with r the area radins and dff^ the standard metric of a nnit 2-sphere. Then the field 
eqnations and the wave eqnation for 0 = 0 (m, r) can be written as 


(In^),^ = r(0,^)^ (7) 

= (1 - 2rV(0)) c/, (8) 

g{g/g),u = 2r - #,«0,r], (9) 

r~^{r‘^g4>^r),r = ‘2(j)^u + 2r0,™ + rgV\(j)). (10) 

For the center r = 0 to be regnlar (finite Ricci scalar), we take 

^(m,0) = ^(m,0) = 1, (11) 


which also fixes the nnll coordinate u as the proper time of an observer at the center. 
Additionally, the Misner-Sharp (Hawking) mass is defined by 


2m{u, r) 


1 - 


g 

g 


( 12 ) 


To study self-similar solutions, we notice that including an exponential potential as given 
in (|2]) does not disrupt the existence of a homothetic Killing vector {L^g^^ = 2gfj_y) 


^ = udu + rdr 


(13) 
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for the metric 
scalar field (see 



This allows us to adopt the same self-similar ansatz for the massless 
and the Appendix A of for detailed derivation): 


ds^ = —g{x)g{x) dn^ — 2g{x) dudr -|- 

r 

(j) = hix) — fi:ln(—u), x = -, u < 0 , 

u 


(14) 

(15) 


with h a real-valued unknown to be determined. Apparently, the self-similar coordinate x 
becomes singular at u = 0 , and we are primarily concerned with the past region u < 0 as u 
increases toward the future. Moreover, note that the constant k in the ansatz also serves as 
a parameter in the scalar potential (E]), and hence is hxed by the Lagrangian if Vq 7 ^ 0. One 
can check that this particular arrangement is necessary in order to reduce the held equations 
to ordinary differential equations (ODEs). 

Derived from dsnQi), the resulting equations for the four unknowns— g{x), g{x), h{x), 
and 7 ( 0 ;) (dehned in flT 6 |l i— consist of four hrst-order ODEs ( ' = d/dx) plus one algebraic 
equation: 


xh' = 7, (16) 

xg' = ^7^ (17) 

xg' + g= {I-2V,x\^~^/^)g, (18) 

{g — 2x) (xy) = 2kx — {g — 2x) 7 -I- 2VQX^y^l*^gy -|- (19) 

g — g = 2fyx — {g — 2a;) (7^ -I- 2^7) -f- 2 Vi^x^y^l'^g. (20) 


One may check that the derivative of the last equation, which derives from (E]), is a combi¬ 
nation of (llOtllOp . If Vo = 0, the equations simplify to the special case for massless scalar 
helds (see 0, (2.9-2.13)). 

To convert the above equations into an autonomous system, we follow Brady’s method 
and introduce new unknown variables: 

w(s) = x^yyi'^ > 0 , yis) = -, z(s) = 3 , a; = e® > 0 . ( 21 ) 

9 9 
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Thus we obtain (' = d/ds = xd/dx] cf. [7|, (2.18-2.21)) 


w = 2 t (;(7 + k)/K, 

y=l-y(l+'t'^)- 2V'o«i, 

4 = z [2 - (1 - 2Vow)y^'] , 


( 22 ) 

(23) 

(24) 

(25) 

(26) 


(1 — 22)7 — 2kz — 7 (y ^ — 22 ) + 2V(,wy ^(7 + a ^), 


(1 — 22)(7 + k)^ — 1 + — (1 — 2VQw)y Vo — ±1 or 0. 


It can be verified that 02311261) imply 0221 ) . Hence using 0261) . we can remove the term VqW 
in 023II25P and reduce the problem to a 3-dimensional dynamical system: 


9 = 2(7 + (t) 1(7 + (t)2 - 7] y, 

2 - - 2(7 + k)^2^ + [(7 + fr)^ + 1 - 2, 


(27) 

(28) 
(29) 



Given these equations, w as determined by fl26]) satisfies fl22]) . Other ways to bring about 
autonomous systems are also possible. For instance, see the Appendix □ for comparison with 
Christodoulou’s treatment and notation. 

The system 0271l29p . as derived from 0221l26p . automatically includes Brady’s special case 
for massless scalars, i.e., 022112 6 p with w = 0 (see also |^, (2.18-2.20)). One can retrieve the 
latter from 02 7112 9 p by imposing the algebraic constraint (j26p with Vqw = 0. 

To facilitate analysis of the autonomous system, we further introduce 


This substitution is intended to bring z = -|-cx) to a finite region (compactification), i.e., 
C = 1, as we will see that naked singularity solutions tend to z = -|-cx) as s —)■ -|-cx). In terms 
of y, C, and 7 , the system 02711291) takes the form 


2/(7 k)(37C -2'y + O 


(31) 

(32) 

(33) 



C = - [ 2(7 + k)2 + 1 - K^] + [(7 + Kf + 1- K^] C, 

/ ,2.2a \ 1 A 


and the constraint for massless scalar fields (([26]) with VqW = 0 ) reads 



( 34 ) 
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Furthermore, the initial conditions imposed by regularity at the center (cf. flTT]) l require 


w —)• 0, 1/ —t 1, C 0, 7^0, as s ^ —oo {x 0). (35) 

The objective is to solve the evolutionary equations fl3HI33p from s = —oo (u < 0, r = 0) to 
possibly s = +oo (u = 0,r > 0). Among solutions satisfying fl35|) . black hole formation is 
indicated by integral curves reaching the plane y = 0, where it locates an apparent horizon 
(cf. fll2p h Hence for naked singularities, we need to identify solutions that can avoid 
contacting y = 0 for the entire range of s. For this purpose, solutions terminating at a hnite 
critical or singular point as s —)■ +cxd are reasonable candidates. Moreover, requiring that 
the mass function be non-negative implies y < 1, which can be used to narrow the search. 

In addition, the system fl3HI33p is invariant under the mapping n —)■ —k, 7 —7 with y 

and ( unchanged. Without loss of generality, we only consider k > 0 henceforward. 


III. LOCAL BEHAVIORS 

Before we investigate global solutions, it is important to understand local phase space near 
critical and singular points (see Tabled] for a quick summary). At these key locations, we can 
determine properties such as stability, dimensions of stable/unstable manifolds, uniqueness 
of solutions, etc., to help us build a fuller picture. Because the phase space for the massless 
scalar held is contained as a surface (signihed by the constraint fIMP ) in our 3-dimensional 
phase space, we expect that many discussions on this surface can be transferable to the 
whole space. 


A. The initial point s = —00 
The initial point 

0: y = l, C = 0, 7 = 0, (36) 

is also a critical point of the vector held dehned by fl3HI33p . Therefore, we can determine the 
behavior of solutions in a neighborhood of O by the standard linearization method. With 

y = l + xi, C = X 2 , l = xz, (37) 
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the Taylor series expansion in xi^2,3 gives rise to a linearized version of fl311l33p : 


x = Ax, X = (a;i,a;2,a;3)'^, 


The Jacobian matrix A has eigenvalues 


A 


^ 0 -2k^ 

0 1 0 

0 1 y 


Ai — —1, A2 — 1, A3 — 2. 


(38) 


(39) 


Thus (9 is a saddle point. The general solution of 


xi = cie^i^ 


cse 


A 3 S 


X2 = C2e^2*, 


^3 = 


is given by 

Ole'll® + C2K^e^^^ + 2c3e^®^ 


2k 


(40) 


with three integration constants ci,2,3- Now we recall the conditions (|TT|l . (ET]), and ([30]). 
They together imply 

lim C e“® = 2, (41) 

s —>•—00 

and hence set C2 = 2. Consequently, solutions subject to the conditions fl3^ and flTTl) take 
the form 


y = 1 — + 0{xiXj), ( = 2e^ + 0{xiXj), 7 = /te® + —e^® + 0{xiXj), (42) 

K 

which has one free parameter C3 associated with A3 and i,j = 1,2,3. Particularly, the 
massless solution, which satishes fl3T|) . occurs uniquely at C3 = 7k^/3. The presence of a 
free parameter can also be conhrmed by the Taylor series solution at O in terms of x. More 
specihcally, we have 


y = 1 — ctx^ — 2ctx^ + 0{x^), 


( = 2x — 4a:^ + (2ct + 8 — k‘^)x^ + 0{x'^), 


(43) 


'y = Kx + [ — + k]x^ + 


3ct / ,, 

—— + (ct + 7)k — — 

Kj Zi 


x^ + 0(a:^), 


with Ct = k^/ 3 for the unique massless solution. These series solutions make up our hrst 
step toward global solutions. They will also be useful for generating approximate initial 
values near O for numerical integration. 

Compared to Christodoulou’s massless case, the interior solution is no longer unique. 
Instead, we have a 1-parameter family of integral curves emitting from O. This enlarged 
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solution manifold is due to the addition of the third eigenvalue A3 (Ai_2 identical with the 
massless case j^), which will lead to richer behavior. 

So far, we have not specihed the sign of the scalar potential since the factor Vq does not 
appear in (13111331) . To determine the sign, we plug the series solution fH3|) into (1^ and 
obtain 

0 < = ^ = + 2x) + 0{x^). (44) 

/I/q 

For a real scalar field, the above quantity must be kept non-negative, which then sets 

f —1 if ct < k^/3 (cs < 7k^/3), 

1/0 = <^ ^ ^ ’ (45) 

y +1 if Ct > fi:^/3 (cs > 7 k^/3). 

Likewise, using (j42]), we have an equivalent sign choice in terms of C3 given in the parentheses 
above. Therefore, based on the sign of Vq, the 2-dimensional interior solution manifold is 
divided into two regions by the unique massless solution (cr = «;^/3, Vq = 0) (cf. Fig. [T]). 


B. Critical points s = -|-oo 


Besides the initial point, the system possesses three more critical points that may serve 
as ending points of solution curves. We can treat them by linearization as well. 

For the first critical point 


Vi : y 



2 

, 7 

o ~\~ Kj 


1, 


the matrix A reads 


with eigenvalues 


A 


0 ^(3 + k)^ 

0 -2(1+ k) 

4(1+k) (1+k+k^)(k+3)^ 

\^k(1—k) k(1—k) 


— K 


\ 


4(1-k)(1+k) 
(3-hK)^ 
2(1+^— 
k(1—k) J 


K — \/4 — 3k^ 

1 — K ’ 


A2 


K A \/4 — 

1 — K 


A3 


2(1+ k) 


(46) 


(47) 


(48) 


Particularly, we have ReAi^2 < 0 and A3 > 0 when k > 1, and Ai < 0, A2 > 0 and A3 > 0 
when 0 < K < 1. Thus Vi is a saddle point. For the second critical point 

'P-i- y = \^ C = 7 = -1, (49) 
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the matrix A reads 


A = 


0 

0 

\k(1+k) 


|(3 — —k ^ 

-2(1-«) 

(1—k+k^)(k—3)^ 2(1—K+K^) 

/^(i+k) k(i+k) y 


(50) 


with eigenvalues 

Ai = - 


-K 


\/4 — 3k^ 


Ao — 


-K 


+ V4 - 3fi;2 


A,-? — 


2(1 - k) 

K 


(51) 


1 + K 1 + K 

Particularly, we have ReAi,2 < 0 and A3 > 0 when k > 1, and Ai < 0, A2 > 0 and A3 < 0 
when 0 < K < 1. Thus V-i is also a saddle point. 

Both of these points reside in the massless phase subspace, and they are directly inherited 
from Christodoulou’s case for having the same locations as well as the first two eigenvalues 

n 

Ai^ 2- Hence we have denoted them in the same way as in 

The third critical point does not have a massless counterpart: 

1 


Vr 


y = 


C = 0, 7 = -K. 


(52) 


1 — 

However, as we will see in later sections, it is not of much relevance to our discussion. Here 
we only list the matrix A and its eigenvalues below without further ado: 


/ 


A = 


V 


0 0 

0 l-K^ 


2k 

1—k'^ 

0 




0 


K 


Ai,2 — 


T — =F a/(1 + — 7) 


A3 — 1 — fi) 


(53) 


(54) 


C. Singular points C = 1 ? •? = +00 

Besides the points V±i , Christodoulou also discussed a third critical point 

^0 : y= . ^ 2 , C = 1, 7 = -«, (55) 

1 ~r 

in the 2-dimensional phase space, which turned out to be very important since solutions 
representing naked singularities tend to this point as s ^ -f-cxo. For our 3-dimensional 
problem, Vg expands into a straight line in the ^/-direction: 

'Po(yo) ■ y = yo^o, c = i, 7 = -«• (56) 
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Points on the line are critical for the equations fl321l33p . but singular for (1^ . 
To probe into the local behavior near Vo{yo), we continue using 


y = yQ + xi, ( = l + a;2, 7 = -k + X3, 


(57) 


as in previous subsections and perform (Laurent/power) series expansions in terms of Xi^2,3- 
Only keeping terms up to the hrst order, we arrive at 


oc 

Xi = 2yoKX3 - I/O—, 

X2 

X2 = - 1)X2, 

X3 =-h X 2 - 3:3. 

I/OK / 


(58) 

(59) 

(60) 


It is noted that a partial linearization succeeds in the {(, 7)-plane since the linear equations 
and fl60|l do not depend on xi. The matrix A for X2,3 alone reads 


K 


2-1 0 


+ -1 

yoK 


(61) 


The eigenvalues are 


Ai — —1, A2 — —(1 — ^2), 


(62) 


which coincide with the massless case (cf. j^, (3.16)). Solving the linear equations gives rise 
to 


C = 1 + C2e H-, 7 = -K + cie * + C2 


1 - 1/0 


+ ^] e-P-^b + ... , (63) 


with two parameters ci^2- Fortunately, the troublesome term x|/x2 in 
under 0631) . and we obtain formally 


is well-behaved 


, C2(1-|/o + |/ok^)(1-|/o-|/ok^) -(i-k2)s , 2ci(l-|/o) 
y = yo + -^7^-77-e ^ ^ ^-e 

K°(l — K )|/o K 

I yo4 (iw)s I 

C2(1 + «)2) ^ ’ 


(64) 


where quadratic and higher order terms in Xi^2,3 are omitted. Hence for < 1, the asymp¬ 
totic expansions suggest that each Po(l/o) acts like an attractor (except in the //-direction). 
This will be later conhrmed by numerical calculation. Very much like the massless solutions, 
'Poiyo) will serve as the endpoints of naked singularity solutions (see Sec. IV)) . 
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D. Singular (connecting) points C = 5 ? s = 


If a continuous integral curve from O were to reach a point Vo{yo) with (^ = 1, it would 
have to cross the plane ( = 1/2 at the following set of singular points: 

y* 


Vsiy*) : y = y^y^0, C = 7 = - . 

2 


(65) 


which is imposed by fl3^ (multiplied by 1 — 2/' on both sides). Hence if an interior solution 
attains Vs{y*) at some hnite s = s*, it is necessary to extend it past C = 1/2. 

To tackle these singular points where the uniqueness of solutions fails, we follow Christodoulou’s 
method and adopt a new independent variable t such that 

ds 


dt 


= -2(l-2C). 


( 66 ) 


In terms of t, the singular point Vs{y*) becomes a critical point of the transformed equations. 
Then perform linearization with 


y = I/* + xi(t), 

and calculate the matrix A as 

^ 0 
0 


C=2+^2(t), 7 = 


y*Kp 

l -«2 


y* 






(67) 




v? 

along with its eigenvalues 


1 

. .7 Q 


( 68 ) 


Ai — A2 — 1 — A3 — 0, 


(69) 


which are independent of y*. The zero eigenvalue A3 reflects the fact that the set of VAy^) 
forms a curve, while the other two eigenvalues are identical to the massless case (cf. |5|, p. 
620). Therefore, if < 1, both Ai^2 are positive and one can construct a 1-parameter fam^ 
of solutions which tend to Vs{y*) as t —)■ —00 (or s —)■ s* from either side of C = 1/2) |5[. 
Consequently, a single interior solution curve from ( < 1/2 may branch out at Vs{y*) into 
a 2-dimensional solution manifold beyond s*. We will further test this claim by numerical 
calculation. 

As an extra comment, we mention that the massless phase subspace determined by 
only intersect the curve Vs{y) at one point: 


Vs 


y = 


1 + 


C = o’ 7 = -• 

2 K 


(70) 
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which was closely examined by Christodoulou j^. In Sec. IIVBI we will compare continua¬ 
tions of the solutions from Vsc with those from other points on Vs{y)- 

E. Apparent horizon y = 0 


The apparent horizon is signaled hj y = 0, which also renders the equation (133|1 sing ular. 


To capture the dominant behavior near y = 0, we perform the Painleve analysis 
identify the following Puiseux series expansions: 

/4\/2aiK ^ \/2(2aia2 — 1) 


29 1 and 


y \/®oo 5 


— -1“ 4aia2 ± . 

002 \ 


3a2K 


C — 2 ®2\/Soo ~ S + • • • , 


7 = ± — 

v2\/’Soo - s 

which converges to 


2k, 20x02 — 1 


\/2k^ 


601O2K 12 o2 


+ 


Vah = P.IO) : 9 = 0, C = 7 = ±70, 


(^co — 5 ) + ■ ■ ■ , (71) 

(72) 

(73) 

(74) 


as s —)■ Soo—. This series solution is general for having three free parameters Oi, 02 and Soo- 
Thereby, one may regard the apparent horizon Vah effectively as an attractor in the phase 


space. The leading order terms of the expansions agree with Christodoulou’s Proposition 
3.2 in js]. 

Using d73|) . we can estimate the scalar field near the apparent horizon as 




x\/\n xah 


dx =-4 ^In^^^^W (75) 

— In X V V X / 


with \axAH = Soo and x < xah- Hence, the scalar field maintains finite and continuous at 
the apparent horizon and may be further extended to infinity. 


IV. GLOBAL SOLUTIONS 

We move on to constructing global solutions via numerical integration and verifying our 
previous analysis on local behaviors. To systematically create numerical curves from the 
initial (stationary) point C>, our strategy is to use the truncated series solution flT2]) or 
fl43|) to generate initial values near O (e.g., taking s = —10 and C3 = 1,2, • • •), and then 
integrate both backward and forward with the standard fourth-fifth order Runge-Kutta 
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TABLE I. Summary of important points in the phase space. 


Points Coordinates (y, C, 7) Description Stability Eigenvalues 


0 

(1,0,0) 

initial starting point s = —00 (tt < 0, r = 0) 

saddle 

(IMl) 

Vi 

^ 1) 

V2’ 3+k’ 

critical ending point s = -|-cx) (u = 0, r > 0) 

saddle 

iS]) 

V-I 

(1 ^_ 1) 

O’ 3-k’ D 

critical ending point s = -|-oo 

saddle 

m 

Fo(y) 

(y + 0 , 1 ,-k) 

singular ending points s = -|-oo, k G (0,1) 

attracting 

(IMD 

T^siy) 


singular connecting points s = s*, k € (0,1) 


(EH) 

'Psc 

( 1 1 1\ 
''T+iF’ 2’ K> 

singular connecting points s = s*, k G (0,1) 


(EH 

Vah 

(0, ^,± 00 ) 

apparent horizon s = Soo 

attracting 



method. This simplihes numerical procedures and fulfills the conditions fl35ll and flTTl) to a 
satisfactory accuracy. Moreover, upon nearing the singular curve Vsiy) or the plane y = 0, 
it is also desirable that step sizes of integration be kept small so as to produce more accurate 
plotting. Similar procedures also apply to integral curves from other critical points. 

As we have seen in local studies, the nature of the critical and singular points relies 
heavily on the value of the parameter k. Our numerical results split into two cases. 


A. Case k > 1 

In our first set of plots with k = 1.2, solutions of interest are limited to 0 < C < 1/2 
with Vi involved. Since Vi is a saddle point, the 1-dimensional unstable manifold (red 
curve) can draw away incoming integral curves (green and cyan) from O on either side of 
the 2-dimensional stable manifold (cf. Sec. IIIIBll . This naturally classifies the solutions 
into three types. For solutions heading toward the plane y = 0 (green curves), they have 
ct > and correspond to Vq = 1 by (H5|) . Numerical integration also indicates that they 
quickly exceed the maximal bound of the programme within a finite s while ascending in 
the direction of Vah (cf. (17^ ). This hrst type represents black hole formation. On the 
other side of the stable manifold, solutions moving away from y = 0 (cyan curves) have 
Vo = — 1 and are deemed unphysical in the sense that they severely violate the positive 
mass condition y < 1 at a large s (in fact, the mass function can go arbitrarily negative). 
The intermediate between the two types is an exceptional solution (blue curve) terminating 
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FIG. 1. Numerical integration with k = 1.2. The right plot shows the back view of the left. Green 
{ct > k^/3, Vo = 1) cf. (US])), blue, and cyan (cr < k^/3, Vq = —1) curves emanating from O 
are members of the 1-parameter family of solutions of the system (13111330 subject to (|35p and (|4ip . 
The negatively unstable manifold of O and the unstable manifold of Vi are marked by red curves. 
The blue curve {ct = k^/3, Vq = 0) ending at Vi delineates the unique massless solution. The 
maroon curve represents the set of singular points Vs{y) with C, = 1/2. 


5|]. According to j^, the 


at Pi as s +00. It also lies in the massless phase subspace determined by (BH), hence 
corresponding to the unique interior solution by Christodoulou 
associated spacetime develops a central singularity at m = 0 that is not visible to observers 
at a hnite radius. 


Unlike the massless evolution, black hole formation completely takes over under positive 
exponential potentials. This “regularizing” effect appears in line with cosmic censorship. 
Nevertheless, negative exponential potentials seem to have an opposite effect by rendering 
the spacetime more pathological. This comparison will be further discussed in later sections. 
For other values of k > 1, the qualitative behavior we have described remains the same. 
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FIG. 2. Numerical integration with k = 0.5 (roughly front and back views). Integral curves from O 
may reach the singular curve Vs{y) (maroon). Particularly, the blue (Vq = 0) curve attaining Vsc 
is the unique massless interior solution interpolating between green {Vq = 1) and cyan {Vq = —1) 
curves. In the region C, > 1/2, curves drawn from V± and Vsc are continuations of the interior 
solutions. The two purple curves terminating at V-i separate naked singularity solutions and black 
hole solutions. The red curve passing through Vi iV-i) with one end reaching Vo{y) represents 
the unstable manifold of Vi (P-i) in the massless phase subspace 


B. Case 0 < k < 1 

This case is much more involved as integral curves can pass through the plane ( = 1/2. In 
the region 0 < (/ < 1/2, an interior solution curve shown in Fig. |2]and Fig. Elcan either move 
toward the apparent horizon y = 0 or run into the singular curve Vs{y) in a hnite s. In Sec. 
IIII Cl we point out that once attaining a point on Vs{y), an interior solution may branch into 
a one-parameter family of exterior solutions. This is conhrmed in Fig. [2] and Fig. [3], where we 
have picked out three representative locations on Vs{y) — Vsc for massless solutions and V± 
for Vq = ±1—from which continuations are plotted. Amongst these extended curves, naked 
singularity solutions correspond to those approaching Vo{y) as s ^ -foo. Clearly in Fig. [2] 
with K = 0.5, such curves can emerge from Vsc (also V-, but not from V+). However, when 
K = 0.7 as in Fig. [3l they cease to show up in massless solutions altogether (blue curves). 
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FIG. 3. Numerical integration with k = 0.7 (roughly upper front and back views). Behavior of the 
green (Vq = 1) and cyan (Vq = —1) curves are similar to those in Fig. [2l However, for massless 
solutions, the blue curves and the unstable manifold of Vi (red curve) no longer extend to the 
singular line Vo{y), indicating that naked singularities cease to occur. 


This numerical result agrees with Christodoulou’s rigorous argument that, in massless scalar 
fields, naked singularities only exist when < 1/3. To understand this phenomenon more 
intuitively from our plots, we note that the unstable manifold A4i (red curve, restricted on 
the surface fIM)) ) of the saddle point Pi appears to be a good indicator since it may dictate 
asymptotical behavior of the integral curves in its vicinity. More specifically in Fig. [21 one 
arm of extends to Vo{y), thereby pulling solution curves nearby to the same endpoint. 
While in Fig. |21 both ends of A4.i lead away from Voiy), and are followed by other integral 
curves. Moreover, according to numericsMlie critical case of k = 1/ y/S has the lower segment 
of Ml falling exactly on the point P_i No naked singularity arises in this critical case. 

The situation, however, differs greatly for solutions with negative exponential potentials 
(cyan curves). Starting from a singular connecting point lower than Psc, an exterior solution 
may evade the influence of Pi and still reach out to Po(y) regardless of < 1/3 or 1/3 < 

< 1. The only caveat here is that solutions of such may locally violate the positive 
mass condition y < 1 and may even have negative total masses (see Fig. [3]). A similar 


situation was also encountered in 


18j. Besides naked singularities, other outcomes are also 
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presented. In both Fig. [2] and Fig. [3l by tnning the second parameter, one can observe 
a transition between naked singularity and black hole solutions. The boundary separating 
these two extremes is demarcated by critical solutions terminating at P_i (purple curves; 
more generally, the stable manifold of V-\). The resulting spacetime can be thought as 
sitting on the threshold of containing a naked singularity or a black hole, but just avoiding 
both. Similar critical behavior is also found in massless solutions (blue curves near Vi and 


V-i in Fig. [2]). As an aside, see 30, iMj for more examples on naked singularity/black hole 


transitions with non-trivial scalar potentials. 

As regards solutions with Vq = 1, black hole formation appears to be the only possible 
outcome as all green curves are heading to Vah instead of Vo{y). This observation is further 
supported by numerical tests on other locations on Vs{y) above Vsc- Therefore, we speculate 
that the naked singularity solutions by Christodoulou defy extensions to positive exponential 
potentials. 


V. SPACETIME STRUCTURE 

In this section, we examine various aspects of the metric (1T4|1 with inputs from the 
solutions we have computed. To begin with, we show that for solutions terminating at 
Vo{y), the corresponding spacetime indeed contains a curvature singularity at the point 
n = 0, r = 0 that is visible to distant observers. To see this, one can calculate the Ricci 
scalar as 


R = -2T = 2(V„0V“0 + 4U(0)) 


= 2 


— + -(0,r)^ + 4Uoexp 

9 9 


2g 


u^x^g 


2?/ 


7 


9 
2x , 

9 

c 

1-C 


7 ^ - — ( 7 ^ + « 7 ) 
9 


+ 


(?) 

8 Vote 




(7^ + «7) 


+ 


SVqW 

u^x'^ 


(76) 

(77) 

(78) 

(79) 


and the limit y ^ y^ as s ^ +cx) (equivalently 


Then by applying the asymptotic series 
u —)• 0—), one arrives at 

R=^ 2{yo-l)(^^-^-Ayo^ , (80) 

with I/O = 1/(1 + for fhe massless solutions. Clearly, this scalar invariant blows up at 
r = 0, but remains hnite elsewhere on the future light cone m = 0 of the center (the Cauchy 
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horizon). Therefore, the spacetime may be further extended beyond m = 0, at which the 
original self-similar coordinate fails. 

Now it remains to establish that the central curvature singularity is indeed visible, i.e., 
that the outgoing null geodesics u = 0 can escape from r = 0. Following the method in 
(see also js]), we investigate incoming radial null geodesics, assuming ( > 1/2, and show 
that they can approach a hnite r as u —?• 0. In terms of the variables t = — ln(—u) and 
s = ln(—r/u), the null geodesic equation 


dr 

du 


I 

'2 


can be rewritten as 


ds 

dt 


Cis) 


2 > 0 . 


Hence, along a null ray originated at some (to, so > s*), we have 

c 


The asymptotic series 


t-to = 
implies that 


-50 


2C-1 


ds'. 


(81) 


(82) 


(83) 




s —)■ -|-CX). 


(84) 


Therefore, the quantity 


In r = s — t = So — to — 


'so 


c 


2C-1 


- 1 ds' 


(85) 


converges to a hnite limit as t —)■ -|-cx) since the integral above is bounded. Also because 
C2 < 0 from numerical results, we conclude that along an incoming null ray near u = 0, the 
radius r decreases to a hnite value as u increases to 0. 

Using the limit y ^ yo as s ^ -|-cx), we can further determine the mass function at u = 0 
as 


m = 


yo 


( 86 ) 


with yo = 1/(1 -|- for the massless solutions and yo > 1/(1 + «(^) for Vq = —1. Thus the 
central singularity itself has vanishing mass. Similarly, we can write down the scalar held in 
the limit u —)■ 0 as 


= 2^" 


1 - (1 -F K^)yc 

2Vo 


- nlnr, Vq =-1, yo>l/{l + K^), 


(87) 
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and for the massless scalar fields, 


(j) = —K Inr, (88) 

both of which diverge logarithmically at r = 0. From the equation flHTD . the negative 
potential energy at n = 0 obeys 


^oc--, 


(89) 


which is non-negligible near the center. 

Lastly, we comment that the naked singularity solutions that we have considered are, by 
themselves, neither asymptotically flat nor AdS. Such unusual asymptotic is in fact typical 
for negative exponential potentials 
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34j |. Nevertheless, we should emphasize that since the 


past light cone of the singularity is given hj —r/u = exp(s*), the naked singularity formation 
is essentially a local event near the symmetry center and independent of the geometry far 
from it. 

VI. CONCLUDING REMARKS 


In the class of self-similar spacetimes, we have studied the effect of exponential poten¬ 
tials on the naked singularity solutions that were hrst discovered by Christodoulou and 
numerically explored by Brady. For every hxed k, the original 2-dimensional phase space is 
embedded into a 3-dimensional one, which allows more dynamical events to take place. Us¬ 
ing numerical methods, we build a pictorial representation of Christoloulou’s solutions and 
put them in a richer context with new solutions. In our 3-dimensional plots, two competing 
“attractors”—one for naked singularities and one for black holes—play a major role as a time 
evolution may be drawn into either one of them. On one side of Christodoulou’s solution 
manifold, a two-parameter family of naked singularity solutions continues to exist under the 
negative potentials V = — exp{2(f)/ k) for 0 < < 1. In contrast, on the side with “reg¬ 

ularizing” positive exponential potentials, all numerical evolutions tested by us encounter 
apparent horizons, which suggests that the naked singularity solutions by Christodoulou 
may be unstable against perturbation of the positive potentials. A similar effect is also 
found in the regime > 1 (see Sec. IIV A|l . 

Regarding future work, it would be natural to investigate embedding of our solutions in 
higher dimensions 3^ [3^. In a follow-up research 3^, we will show that such embedding 
leads to naked singularity solutions in higher dimensional vacuum, thereby putting our 
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current studies in a fruitful perspective. To conclude, our results may shed light on the issue 
of cosmic censorship in string theory. 


Appendix: Christodoulou’s notation 

One can translate Brady’s notation that we have been using to Christodoulou’s through 

z~^ =2{1 —13), J = 0, n = k. (A.l) 


In terms of the latter, the system (131113311 reads 


ds 

ds 

ds 


{9 + k){2p9 - 9 + k) 

2 ( 1 -/ 9 ) ’ 

1-k^ - [{9+ k)^+ 1- k^] 13, 

^{k9 - 1) + [(0 + k)^ - (1 + P)] 9 + ^ 


l + k^ + 


1-/3 


{9 + k)^ 



(A.2) 

(A.3) 

(A.4) 


Also the constraint (IHTll for massless scalar helds becomes 

e2^ = l + fc2 + _^(0 + fc)2. (A.5) 

Using this equation, we can remove A in (1A.4I1 and obtain a reduced system 

^ = 1 - k'^ - [{9 + k)'^ + 1 - f3, (A.6) 

ds 

39 k 

^ ^ ~ ^) + [(^ + ~ (1 + 

which recover the equations (0.27a,b) in j^. It is straightforward to verify that if the above 
system holds for (3 and 9, X as determined by flA.Sp automatically satishes (IA.2p . 
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